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Abstract 



We discuss properties which must be satisfied by a genetic network in order 
for it to allow differentiation. 

These conditions are expressed as follows in mathematical terms. Let F be 
a differentiable mapping from a finite dimensional real vector space to itself. 
The signs of the entries of the Jacobian matrix of F at a given point a define 
an interaction graph, i.e. a finite oriented finite graph G{a) where each edge 
is equipped with a sign. Rene Thomas conjectured twenty years ago that, if F 
has at least two non degenerate zeroes, there exists a such that G{a) contains a 
positive circuit. Different authors proved this in special cases, and we give here 
a general proof of the conjecture. In particular, wc get this way a necessary 
condition for genetic networks to lead to multistationarity, and therefore to 
differentiation. 

We use for our proof the mathematical literature on global univalence, and 
we show how to derive from it several variants of Thomas' rule, some of which 
had been anticipated by Kaufman and Thomas. 
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When studying interactions in a system of biochemical compounds, it is quite 
rare that one obtains quantitative results. One can show that, in a given tissue, 
(the product of) a gene A is an activator (or a repressor) of the expression of a 
gene B, but the strength of this interaction, the concentrations and their kinetics 
are usually unknown. The resulting information is essentially summarized by an 
interaction graph G, by which one means a finite oriented graph together with 
a sign for every edge. The vertices correspond to the members of the network, 
and there is a positive (resp. negative) edge from A to B when A activates 
(resp. represses) the synthesis of B. Note that there can exist both a positive 
and a negative edge from A to -B, since A can activate B at some concentration 
and repress it at some other. 

These interaction graphs can be quite complicated. It is therefore very desir- 
able to find ways to use them to restrict the possible behaviour of the network 
they represent. In this paper, we shall address the question of when the net- 
work is susceptible to have several stationary states. A beautiful conjecture of 
R. Thomas [13] asserts that a necessary condition for multistationarity is that 
G has an (oriented) circuit C which is positive, i.e. such that the product of 
the signs of the edges of C is positive, at least in part of phase space. This was 
already proved in several cases [11] [12] [7] [3], and we shall present a proof in 
the general case (Theorem 1). 

To formulate the question in mathematical terms, we choose a continuous 
model. Let n > 1 be an integer, and let F be a differentiable map from R" 
to itself. The evolution of a network of n compounds can be modelled by the 
differential system 

where a; is a differential path in M": the components of x{t) are the different con- 
centrations at time t. To ask whether this system has several stationary states 
amounts to ask whether F has several zeroes. Given any a S M", the interaction 
graph G{a) is defined from the signs of the partial derivatives {dfi/dxj){a) of 
the components of = (fi) at a (see 2.1 below for a precise definition) : a 
positive (resp. negative) sign indicates that j is an activator (resp. a repressor) 
of i. The precise formulation of Thomas' rule is then the following: if F has 
at least two (non-degenerate) zeroes in M", there exists a G ffi" such that G{a) 
contains a positive circuit. 

The main remark leading to a proof of this assertion is the following. If F has 
several zeroes, it cannot be univalent (i.e. one to one). Therefore, if we know 
sufficient conditions for F to be univalent, their negation will give necessary 
conditions for F to have several zeroes. And we might then deduce properties 
of G{a) from these necessary conditions. 

It turns out that finding sufficient conditions for F to be univalent is a 
classical issue in mathematical economy, when one wants to know that the 
factor prices are uniquely determined by the prices for goods. Several results 



3 



have been obtained with this application in mind. G. Gale and H. Nikaido [5] 
gave an elegant criterion of univalence in terms of the Jacobian matrix J{x) = 
{dfi/dxj){x): assume that, for every x G M", the principal minors (rcsp. the 
determinant) of J{x) are nonnegative (resp. is positive); then F is univalent. 
As we shall see, this result leads precisely to a proof of Thomas' conjecture 
(Theorem 1). 

Following the same line of arguments, we can also apply variants of the theo- 
rem of Gale-Nikaido, which are discussed for instance in the book of T. Partha- 
sarathy [10]. A motivation for doing so is that one may want to restrict the 
domain of F to be the closed positive quadrant (since a concentration cannot 
be negative). In general, F will be defined on a product fl of n intervals, open 
or not. However, in this generality, there exist counter-examples to the condi- 
tion of Thomas (see 3.5). But, when F has at least two zeroes, we can still 
find properties of its interaction graphs (Theorem 5), as a consequence of a uni- 
valence theorem of Gale-Nikaido ([5] , [10]). Furthermore, it is often the case 
that, when defining the interaction graph of several biochemical compounds, the 
degradation of each of them is not taken into account. We show in Theorem 6 
that this makes Thomas' rule valid in full generality. 

Another refinement of the theorem of Gale-Nikaido, due to C.B. Garcia and 
W.I. Zangwill [6], applies to the case where f2 is closed and bounded. This 
leads to Theorem 7 which, in turn, for general fl, gives some information on the 
location of the zeroes of F (Theorem 8). 

The economist P.A. Samuelson imagined a stronger univalence criterion that 
those above. It is not true for general F, but L.A. Campbell proved it when 
Q, = R" and F is algebraic, e.g. when each component of F is the quotient of two 
real polynomials. In Theorem 9 we translate his result in terms of properties of 
the interaction graph. 

One can expect more graphical requirements for multistationarity. One of 
them is a conjecture of M. Kaufman and R. Thomas [15], which is stronger than 
the original Thomas' conjecture. We have been unable to prove this assertion, 
but we obtained some evidence for it in Theorems 3, 4 and 5. 

The paper is organized as follows. The first section gives definitions about 
graphs, matrices and determinants. It shows basic lemmas, which are standard 
knowledge in the literature on interaction graphs. In section two we define the 
interaction graphs G{a) and we state the conjectures of Thomas and Kaufman- 
Thomas. In the third paragraph we prove the conjecture of Thomas and we give 
some (counter) examples. In section four we discuss the conjecture of Kaufman- 
Thomas. Next, wc give results when the domain fl is not necessarily open, and 
we refine them in section six. Finally, we discuss the case of an algebraic map 
in the last section. 

The literature on the question studied here is scattered in several different 
journals, and it has been crucial for me to receive papers from different authors. 
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1 Graphs and matrices 
1.1 

An interaction graph G = {V, E, sgn) is a finite oriented graph {V, E) together 
with a sign map sgn : E {±1}- In other words, V (the vertices) and E (the 
edges) are two finite sets and each edge e e E has an origin o(e) € V and an 
endpoint t{e) e V (it may happen that o(e) = t{e)). 

A circuit in G is a sequence of edges ei, . . . , such that o(ei+i) = t{ei) for 
all i = 1, . . . , /c — 1 and t{ek) = o(ei). 

A hooping is a collection C = {Ci, . . . , C^} of circuits such that, for all i j, 
Ci and Cj do not have a common vertex. A circuit is thus a special case of 

fc 

hooping. We let V{C) = JJ ^(^0 be the (unordered) set of vertices of C. This 

set V{C) will also be called the support of C. A hooping is called Hamiltonian 
when its set of vertices is maximal, i.e. V{C) = V. Note that hoopings are 
called "generalized circuits", or "g-circuits" in [4], and Hamiltonian ones are 
called "full circuits" in [15]. 

The sign of a circuit C is 

sgn(C7)= []sgn(e)e{±l}. 

When sgn(C) = +1 (resp. —1) we say that C is positive (resp. negative). The 
sign of a hooping C is 

sgniC) = {-ir+\ (1) 

where p is the number of positive circuits in C (cf. [4]). 

Given any subset I C V we let t/ G the interaction graph obtained from 
G by changing the sign of every edge e <E E such that t{e) G /. Given any 
permutation a £ Aut(y) of the vertices, we let aG be the interaction graph 
obtained from G by replacing each edge j — > i by an edge j — > a{i), with the 
same sign. 
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1.2 

Let n > 1 be an integer and A = {aij) an n by n real matrix. We can attach 
to A an interaction graph G as foUows. The set of vertices of G is {1, . . . ,n}. 
There is an edge e with o(e) = j and t{e) = i if and only if a^- ^ 0. The sign of 
e is the sign of a^. 

Given any subset Id {1, . . . , n}, the principal minor of A with support / is 
the real number det (Aj), where Ai is the square matrix {aij)ij^i. By definition 

det(A,)= ^ e((7)nai,(,), (2) 

where S/ is the group of permutations of / and e{a) is the signature of a 
(Recall that e is defined by the equalities e(<T(T') = e{a)e{a') for all a, a' € S/ 
and s{a) = — 1 when a € S/ is a transposition). 
For any a e Sj we let 

a{(j) = e(cr) Y[ ata(i) (3) 

so that 

det (Ai) = ^ a(cr) . 

When a{a) ^ 0, we let sgn(a(CT)) — ±1 be its sign. 

Let D = diag(rfi, . . . , be a diagonal n by n real matrix and / C S any 
subset. It follows from the definition (2) that 

det((^ + £»)/) = ^det(Aj) l[ di. (4) 
Jd iei-J 



Given I c V any subset, we let t/ A be the matrix obtained from A by 

replacing aij by — whenever i G /. Given any a G E„ = Aut(y) we let aA 
be the product of A with the permutation matrix defined by a. Clearly 

G{ti A) = Ti G{A) 

and 

G{aA) = aG{A) . 



1.3 

We keep the notation of the preceeding paragraph. Note that, given any per- 
mutation (T e S/, there is a unique decomposition 

/ = /i n /2 n . . . n 4 
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of / into a disjoint union of nonempty subsets such that the restriction Oa of a 
to la is a cyclic permutation for all o; = 1, . . . , fc. Let C{aa) be the circuit of G 
with edges (i,(Tct(j))i "i & la (note that aio-^(i) ^ since a{<T) ^ 0). We denote 
by C{a) the hooping of G which is the disjoint union of the circuits Gipa)^ 
a = 1, . . . , fc. 

When X is a finite set, we let {X) be its cardinality. The following lemma 
is due to J.Eisenfeld and C.DeLisi ([4], Appendix, Lemma 2), and is probably 
at the origin of the definition (1). 

Lemma 1. Let I C {1, ...,n} he any subset and let a £ T,j be such that 
a{a) 7^ 0. Then the following identity holds 

sgn(C7(a))=sgn(a(<7))(-l)#W+i. 



k 

Proof. Since e(a) = J2 ^{(^a) and / = ]J /„, we get from (3) that 

a=l a 

a{(^) = '[[a{aa) . (5) 

a 

For each a = 1, . . . , A;, since a a is cyclic we have 

The circuit C((Tq,) is positive if and only if Y[ '^iuaH) positive, in which case 
we get from (3) that 

sgn(aK)) = (-l)#(^«)+^ (6) 
When G{(Joi) is negative we get 

sgn(aK)) = (-l)#(^"). (7) 

Therefore, by (5), (6), (7), we get 

sgn (a(a)) = (-1^ (-1) ^ = {-If {-l)*^'^ , 

where p is the number of positive circuits in C{a). Since, by (1), sgn(C(fT)) = 
(— 1)^'+^, the lemma follows. 

1.4 

Prom Lemma 1 we get the following: 

Lemma 2. Let A be an n by n real matrix and let I C {1, ... ,n} be any subset. 
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i) Assume that det{—Ai) is negative (resp. positive). Then there exists a positive 
(reap, negative) hooping in G{A) with support I. 

ii) Assume det(— A/) = 0. Then either there exist two hoopings in G{A) with 
opposite signs and support I, or there is no hooping in G{A) with support I. 

Proof. From (2) we get 

det(-A/) = (-l)#(^)det(A/). 

Therefore, by Lemma 1, 

det{-Ai) = -Y, sgn(C(a))|a(<7)| . (8) 

Assume det(— ^/) is negative (resp. positive). Then, by (8), there exists cr G Sj 
such that C{a) is positive (resp. negative). This proves i). 
If det(— ^/) = 0, either there exist two summands with opposite signs on the 
right hand side of (8), or a{a) is zero for all a G S/. This proves ii), since every 
hooping of G{A) with support I is of the form C{a) for some tr G S/. 

2 Conjectures on multistability 
2.1 

Let fi, be a nonempty interval in IR: 

^i = ]ai,bi[, WiMl, ]aiM] or [aubi], 

n 

with fli > — oo and 6j < +c». Denote by C M" the product =11 ^i- 

i=l 

Consider a map 

which is differentiable, i.e. such that, for each i,j G {!,..., n} and any a G il, 
the i-th component /j of F has a partial derivative ^ (a) at the point a and 




where ||a; — a|| is the norm of x — o and o is the Landau o-symbol (we could also 
use a weaker notion of differentiability, sec [8]). The map F is called when 
every partial derivative dfi/dxj is continuous on f2. 

For any a Gfl, the Jacobian of F at o is the n by n real matrix 



J{a) = J{F){a) = 
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For any a G we let 

G(a) = G(J(a)) 

be the interaction graph of J(a), defined as in 1.2. We also let G{F) be the 
interaction graph defined as follows. Its set of vertices is V = {1,. . . ,n}. Given 
i and j in V, there is at most one positive (resp. negative) edge from j to i; it 
exists if and only if there is a positive (resp. negative) path in G{a) for some 
a e f2. In other words, G{F) is the "superposition" of all the interaction graphs 
G{a). 

Given any subset I C V, we let n F be the map obtained from F = (/j) 

by replacing fi by —fi when i E I. Given tr S I]„, we let aF = {fa—^{i))- The 
interaction graphs of an F are an G{a), a £ Q, and an G{F). 



2.2 

We shall be interested in the set of zeroes of F, i.e. the points n e fi such 
that F{a) = 0. They can be viewed as the stationary states of the system of 
differential equations 

^ = F(.), 

where x{t) is a differentiable mapping from a real interval to f2. 
We say that a zero a of F is nondegenerate when det(J(a)) ^ 0. 



2.3 Conjecture 1. (Thomas [13]) 

Assume that is open and that F has at least two nondegenerate zeroes in fl. 
Then there exists a G ft such that G(a) contains a positive circuit. 

Remark. This conjecture has already been proved in several cases. First, it is 
known to hold when the signs of the entries in ,/(a) are independent of a G 
[11] [7]; see also [12] and the remark in 5.2 below. It was also shown in [3] for 
stable stationary states when O contains the positive quadrant and fi{x) > 
whenever Xi = 0. 



2.4 Conjecture 2. (Kaufman-Thomas [15]) 

Under the same assumption as Conjecture 1, 

i) Either there exist a eCl such that G{a) has a positive Hamiltonian hooping 
and a negative Hamiltonian hooping. 

ii) Or there is a cyclic permutation cr of a subset of {1, . . . , n} and there exist 
a, & G such that the circuit G{a) in G(a) (resp. in G{b)) is positive (resp. 
negative). 

Note that Conjecture 2 for F implies Conjecture 1 for F. It is quite different 
though, since positive and negative circuits play in it a symmetric role. In case 
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i) we shall say that "G has two Hamiltonian hoopings of opposite signs" and 
in case ii) we shall say that "G has an ambiguity". In the latter case, G{F) 
contains two circuits with the same ordered set of vertices and opposite signs 
(but ii) is stronger than that statement). 

3 A proof of Thomas' Conjecture 1 
3.1 

Theorem 1. Assume Q is open and F has at least two nondegenerate zeroes 
in ft. Then, for every I C V and every a £ E„, there exists a £ ft such that 
(jTi G{a) has a positive circuit. In particular, Thomas ' Conjecture 1 holds true. 

Since arj G{x) = G{aTiJ{x)) and since aTi{F) satisfies the hypotheses of 
Theorem 1 if and only if —F does, we just have to check that, for some a G fl, 
the interaction graph G(— J(a)) has a positive circuit. According to Lemma 2 
i), it will be enough to show that, for some a G O, a principal minor of —J{a) 
is negative. In other words: 

Theorem 1'. Assume that, for every a €fl, all the principal minors of —J {a) 
are nonnegative. Then —F can have at most one nondegenerate zero. 

3.2 

For any positive real number A let 

(t)x{x) = -F{x)+Xx. 

Since 

J(</>A) = -J(i^) + diag(A,...,A), 

it follows from (4) that, under the hypotheses of Theorem 1', each principal 
minor of J{4>x) iw positive on il. According to Gale and Nikaido, [5] Theorem 
4, this implies that (px is univalent. The following proposition ends the proof of 
Theorem 1': 

Proposition 1. Let = O ^ R" be a differentiable map defined on an open 
set f2 c M". Assume that, for all A > 0, the map 

<l)x{x) = (j){x) + Ax 

is univalent. Then (j) can have at most one nondegenerate zero in O. 
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3.3 



The proof of Proposition 1 proceeds as [5] Theorem 4' and [10], IV, Theorem 
4, p. 35. Assume a and b are two nondegenerate zeroes of ^ in fl. According to 
[1], XII, § 2.9, p. 477, we can choose open neighborhoods Ua and Ub of a and b 
respectively such that UariUb = 0, Ua C O, Uh C fi, a (resp. 6) is the unique 
zero of (j) in Ua (resp. ?/{,), and the degrees deg((^, a) and deg(0, 0) are 
equal to ±1. Arguing as in [10] (or [5]), loc.cit., we get 

1 = deg(</., C7a U Ufe, 0) = deg(.^, I7„, 0) + deg(<^, 1/6, 0) , 

hence a contradiction. q.e.d. 

3.4 

The restriction to nondegenerate zeroes of F in Theorem 1 is necessary. For 
example, if n = 2, f2 = and 

F{x,y) = {-xy^,-y) 

we have 

Clearly G has no positive circuit. However F{x, 0) = for any a; G M. 
3.5 

It is also essential that be open. Let be the set of {x,y) G such that 
a; > and y > 0. Consider the map F : f2 — > defined by 

F{x, y) = {{y - 2f -x^-l,Ax- 4xy) . 

We get 

Therefore, by (8), G{F) docs not have any positive circuit. On the other hand 
both (0,3) and (0, 1) are nondegenerate zeroes of F in ft. We shall discuss the 
case of an arbitrary O in sections five and six below. 

4 On the conjecture of Kaufman-Thomas 
4.1 

We first assume that n = 2. Let c be as in 2.1 and open, with coordinates 
X and y. For any /i : O — > M we write /i = to mean that h{a) = for every 
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a € fl. If ft. is differentiable, we let hx (resp. hy) be its partial derivative with 
respect to the first (resp. second) variable. 
Let 

be a differentiable mapping. 

Theorem 2. Assume that F has at least two nondegenerate zeroes in fl. Then, 
one of the following conditions holds: 

i) G has two Hamiltonian hoopings of opposite signs; 

ii) G has an ambiguity; 

iii) fxQy = but fx ^0 and gy ^ 0. 

Proof. Assume that G has no ambiguity, and that its Hamiltonian hoopings 
have all the same sign. If, in addition, fxgy ^ 0, since the sign of fx{a) (resp. 
gy{a)) is independent of o, we can multiply / and g by ±1 to get to the case 
where fx{o) > and gy{a) > for all a G f2. It follows that fx{a) 9y{a) > 
and, since all Hamiltonian hoopings have the same sign, fy{a) gxia) < (by 
Lemma 1). Therefore Theorem 1' applies to F and we conclude that F has at 
most one nondegenerate zero. 

Assume now that gy = 0. After multiplying / and by ±1 we can assume 
that /x > and fygx < 0. Once again. Theorem 1' implies that F has at most 
one nondegerate zero, q.e.d. 

4.2 

With an additional assmnption, Conjecture 2 is true for all n > 2: 

Theorem 3. Let F : Q ^ be a differentiable mapping such that f2 is 
open and F has at least two nondegenerate zeroes. Then, one of the following 
conditions holds: 

i) G has two Hamiltonian hoopings with opposite signs; 

ii) G has an ambiguity; 

iii) for any point a e O there exists i G {1, . . . ,n} such that G{a) does not 
contain an edge from i to itself. 

Proof. Assume that all Hamiltonian hoopings in G have the same sign, that 
G has no ambiguity and that, for some a € fi, and for any i & {1, . . . , n}, there 
is an edge in G(a) from i to itself. The last condition means that all diagonal 
entries in J(a) are nonzero. After multiplying each component of F by ±1 we 
can assume that all diagonal entries of J(a) are positive. Let C be any hooping 
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in J(a), and let / be its set of vertices. The disjoint union of C with ah the 
circuits i ^ i, i G {1, . . . ,n} — I, is a. Hamiltonian hooping of G{a). Its sign 

must be the sign of the Hamiltonian hooping which is the disjoint union of all 
the positive circuits i ^ i, i £ {1, . . . ,n}, namely (—1)"+^. We conclude that 

sgn(C) = (-l)#(^). 

By Lemma 1, this implies that C = C{a) with sgn(a((T)) = +1, a G S/. Since 
G has no ambiguity, for any cr € S/ we have sgn (a(cr)) > in fi. Therefore, 
for any x G Cl, all the principal minors of J{x) are nonnegative. Applying 
Theorem 1', we conclude that F has at most one nondegenerate zero. 

4.3 

Concerning the graph G{F) (sec 2.1) wc can prove the following: 

Theorem 4. Let F : fl ^ M" be a differentiahle mapping such that Vl is open 
and F has at least two nondegenerate zeroes. Then G{F) has two Hamiltonian 
hoopings with opposite signs. 

Proof. We first remark that wc can find a G and a E S„ such that none of 
the diagonal entries in J{ijF){a) is zero. Indeed these entries are, by definition, 

— Q^^/" (a), j = 1, . . . , n, and, when a is a non degenerate zero of F, it follows 
from (1) that, for some cr G S„, 

This proves the claim. 

We may then choose / such that all the diagonal entries in J{ti aF){a) are 
negative. In other words, for each i G {1, . . . , n}, the interaction graph t/ aG{a) 
contains a negative edge i ^ i. On the other hand, by Theorem 1, there exists 
6 G r2 such that r/ aG{b) has a positive circuit C. In r/ aG{F) the disjoint union 
of C with the negative edges i ^ i,i (f: V{C), is a positive Hamiltonian hooping, 
when the union of all the negative edges i — > «, z G {1, . . . , n}, is a negative one. 

It is thus enough to show that, if tiG{F) or aG{F) has two Hamiltonian 
hoopings with opposite signs, the same is true for G{F). This is clear for tiG{F) 
and, for aG{F), it follows from Lemma 3 below. 

4.4 

Lemma 3. Let G be any interaction graph, a G Aut{V) a permutation of its 
vertices, and C an Hamiltonian hooping in G. The image ofC in the interaction 
graph aG is then an Hamiltonian hooping with sign e{a)sgn{G). 
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Proof. Recall from 1.1 that aG is obtained from G by replacing each edge 
j — > i by an edge j — > a-(i), with the same sign. As a collection of edges, C has 
a well defined image crC in aG. To check Lemma 3 we may assume that a is 
the transposition of two vertices i and j (transpositions span KutiV)). 

Assume first that i and j are in the same circuit C\ of C. Then all the circuits 
in C other than Ci are fixed by a. The image of Ci consists of two disjoints 
circuits Di and D2- More precisely, if the vertices of Ci are 1 2 ... i ... j ... m 
(as we can assume), we get 

aCi=DiY[D2 

where the sequence of vertices in Di (resp. D2) is 1 2 . . . i — 1 j j + 1 . . . fc 1 (resp. 
+ 1 . . .j — li). Furthermore, if C\ has an even (resp. odd) number of negative 

edges, Di and D2 will have the same (resp. a different) number of edges modulo 
two. From this it follows that aC is a Hamiltonian hooping such that, with the 
definition (1), 

sgn((7C) = -sgn(C), 

as was to be shown. 

Note that we also have 

therefore, by exchanging the roles of C and crC, the previous discussion applies 
also to the case where i and j lie in two different circuits of C. 

q.e.d. 

5 The case of a domain which is not open 
5.1 

We keep the notation of 2.1, where O is an arbitrary product of intervals and 
F : Q ^ R" is differentiable. 

Theorem 5. 

Assume that F is not univalent. Then: 

1) For every cr e S„, 

i) Either there exists a £ fl and i €V such that aG{a) does not contain any 
edge from i to itself; 

ii) Or, for any subset I C V, there exists a G fl such that t/ aG{a) has a 
positive circuit. 

2) When the condition 1) i) above is not satisfied, G{F) contains both a positive 
and a negative Hamiltonian hooping. 

Proof. The map F is univalent if and only if t/ aF is. So, to prove 1), we can 
restrict our attention to —F. Note also that i) is equivalent to the assertion 
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that there exists I (Z V such that aG{a) does not contain any hooping with 
support I. According to Lemma 2 ii), if none of the conclusions in 1) is true, all 
the principal minors of the Jacobian matrix of — are positive on Q. By the 
Gale-Nikaido theorem, [5] Theorem 4, this implies that —F is univalent. This 
proves 1). 

To prove 2), by replacing F hy tj F for an appropriate choice of /, wo can 
assume that, for each vertex i in G{F), there is a negative edge from i to itself. 
Since G{F) also contains a positive circuit by 1), we get, as in 4.3 above, that 
G{F) contains both a positive and a negative Hamiltonian hooping. 

Remark. Since the condition 1) i) is often satisfied, when O is open Theorem 5 
is much weaker than Theorems 1 and 4. 

5.2 

Here is a variant of Theorem 5. 

Theorem 6. Assume given, for every a a diagonal matrix D{a) with pos- 
itive entries. Assume F is not univalent. Then, for some a Gfl, the interaction 
graph 

H{a) = G{J{a) + D{a)) 

has a positive circuit. 

Proof. Assume that, for any a € Cl, none of the circuits of H{a) is positive. 

We know from Lemma 2 i) that all the principal minors of — J{a) — D{a) are 
nonnegative. From (4) this implies that all the principal minors of — J {a) are 
positive, and, again by [5] Theorem 4, F must be univalent. 

Remark. Theorem 6 applies to the situation considered for instance in [14] 
and [9], where —D(a) comes from "terms of decay", which are not taken into 
account when drawing the interaction graph. It was proven by Snoussi [12] 
when the signs of the entries of J(o) are constant. 

6 On the location of stationary states 
6.1 

n 

Assume that $!=]"[ [ai,6i] is a closed bounded subset of R". In that case, 

?=i 

Garcia and Zangwill got a stronger result that Gale-Nikaido ([6], see also [10], 
V, Theorem 1, p. 41). For any I dV define I'^ = V — I and, for any n by n real 
matrix A, let 

mi {A) = det(>l/c) 

and 

mi{A) = nHiy{A) 
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for each i G V. When x lies in fl, we write mi{x) for mi{J{x)) and mi{x) = 
mi{J{x)). Define I{x) C F as the set of vertices i such that Xi = Ui or Xi =bi. 
The result of Garcia-Zangwill is the following 

Theorem 7'. Assume F is and that, for every a £ il, and every subset 
I C lia), 

mi{a) mi(o) > . 

Then F is univalent. 

o n 

In particular, when a £il = Y\]ai,bi[, I{a) = and the only assertion made 

i=l 

is that det J (a) > 0. 

When F is only differentiable, see [8], p. 930, Remark. 
6.2 

Theorem 7' implies the following refinement of Theorem 5: 

Theorem 7. Assume that n is bounded and closed, and that F : CI ^ K" is 
C^. If F is not univalent, for every a € "Sn, one of the following conditions 
holds true: 

i) There exists a € ft and I C I{a) such that no hooping of aG(a) has support 
I'^ , or aG{a) contains two hoopings with support I'^ and opposite signs. 

o 

ii) There exists a G 90 = fi — fi, / C /(a) and hoopings Ci and Ci, for each 

i G I, in aG{a) such that #(/) > 2, the support of Cj (resp. Ci) is I'^ (resp. 
{lY), and 

sgn(C7)nsgn(Q) >0. 

is/ 

Furthermore, either #(/) is even or there exist b G 9f2, /' c I{b), Cji, Cj, 
j & I', with similar properties as above and 

sgn{Cp)l[sgniCj)<0. 
jei' 

Proof. Again, it is enough to treat the case a = 1. When F is not univalent 
we know from Theorem 7' that there is a G f2 and / C I{a) such that 

Assume mi{a) = 0. Then, by Lemma 2 ii), the statement i) must hold. 
Assume now, that, for all J C I (a), mj{a) ^ 0. Since 

m/(a) mi{a) < 
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we must have #(/) > 2. Furthermore, for every vertex k G V,if we replace F by 
TkF = T^k} F the quantity mi{a) is multipUed by +1 (resp. —1) if k G I (resp. 

k e I"). Therefore, in all cases, mi{a) W mi{a) gets multiplied by (— It 

is thus invariant if and only if #(/) is even. If this number is odd, we apply the 
same discussion to —F and we get b and /' such that 

mr{b) Jl mj{b) > 0. 
jei' 

Using (8), the statement i) follows. 
6.3 

n 

Lot il = J]^ rti be an arbitrary product of intervals as in 2.1 and let F : fl ^ K" 

i=l 

be a map. From Theorem 7 one gets some information on where two zeroes 
of F can be: 

Theorem 8. Fix cr e S„. Assume that there exist two points a and b inQ, such 

o 

that F{a) = F{b). Assume furthermore that, when x lies in fl, all hoopings of 

aG{x) are negative. Then there exists x such that, when Xi e dQ.i, the i-th 
coordinate of a orb is equal to Xi, and a subset I C I{x) such that no hooping 
in cfG{x) has support equal to P. 

Proof. One can find a bounded closed product of intervals O' c O containing 
a and b and such that, whenever x & Q.' and Xi e dQ.i, the i-th coordinate of 
a or 6 is equal to a;,. Since the restriction of F to fi' is not univalent, we can 
apply Theorem 7. We are not in case ii) because, J{x) being continuous in x, 
for every a; e O' all the hoopings of G{x) are nonpositive. 

Therefore Theorem 7 i) holds true for some x G Cl', hence the conclusion. 

q.e.d. 

7 The algebraic case 

Assume now that fl = K" and that F is and algebraic, by which we mean that 
its graph {{x, F{x))} c M" xM" is the set of zeroes of a family of real polynomials 
in 2n variables. This will be the case for instance when each component of 
F is the quotient of two polynomials in n variables. In that case, a result of 
Campbell [2] leads to a stronger conclusion than Theorem 5 1). 

Theorem 9. Choose any ordering ofV. Let 



17 



be a map which is not univalent. Then we can choose k < n such that, if I 
consists of the first k vertices, 

i) Either G{F) has two hoopings of opposite signs and support equal to I; 

ii) Or there exists a e K" such that none of the hoopings of G{a) has support 
equal to I . 

Proof. For any k < n and a € R" we let 

dk{a) = det( J(a)/) . 

According to [2], if dk{a) > for all k = l,...,n, the map F is univalent. 
Therefore, under our assumption, dk{a) < for some a and some k. When 
dk{o) = 0, cither i) or ii) is true (by Lemma 2 ii)). Given i € V, when we 
replace F by r^F, d^ia) gets multiplied by —1 (rcsp. +1) ii i < k (resp. 
i > k). So we can assume that, for every I C V , there exist k and a such that 
T/ dk (a) < (with obvious notation) . By replacing F hy n F we see that we 
can also assume that there exist k' and b such that tj dk'{b) > 0. 

Unless i) or ii) holds, we can assume (using Lemma 2 again) that for every 
I CV and every k < n the sign of tj dk{x) does not depend on a; e M". Let 
then m be the maximum of all integers p such that, ii i < k < p, dk (x) and di {x) 
have the same sign. When m < n, we can increase m by m + 1 by replacing F 
by Tm+i F. By repeating this process, we find I such that t/ dk{x) has a fixed 
sign for every k < n and every x € K". As we saw in the previous paragraph, 
this cannot happen. q.e.d. 



References 

[1] p. Alexandroff, H. Hopf, Topologie, Band 1, Berichtigter Reprint, 
Die Grundlehren der mathematischen Wissenschaften, Band 45, Berlin- 
Heidelberg-New York, Springer- Verlag, XIII, (1974), 636 p. 

[2] L.A. Campbell, Rational Samuelson maps are univalent , J. Pure Appl. 
Algebra 92, No.3 (1994), 227-240. 

[3] O. Cinquin, J. Demongeot, Positive and negative feedback: striking a bal- 
ance between necessary antagonists, J. Theoretical Biology 216 (2002), 229- 
241. 

[4] J.Eisenfeld, C. DeLisi, On conditions for qualitative instability of regulatory 
circuits with application to immunological control loops, in Mathematics 
and Computers in Biomedical Applications, J.Eisenfeld and C. DeLisi eds., 
Elsevier, Amsterdam (1985), 39-53. 



18 



[5] D. Gale, H. Nikaido, The Jacobian matrix and global univalence of map- 
pings, Math. Annalen 159 (1965), 81-93. 

[6] C.B. Garcia, W.I. Zangwill, On univalence and P-matrices, Linear Algebra 
Appl. 24, (1979), 239-250. 

[7] J.-L. Gouzc, Positive and negative circuits in dynamical systems, J. Biol. 
Sys. 6 (1998), 11-15. 

[8] M.S. Gowda, G. Ravindran, Algebraic univalence theorems for nonsmooth 
functions, J. Math. Anal. Appl. 252, No.2, (2000), 917-935. 

[9] M. Kaufman, R. Thomas, Model analysis of the the bases of Multistation- 
arity in the Humoral Immune Response, J. Theoretical Biology 129 (1987), 
141-162. 

[10] T. Parthasarathy, On global univalence theorems, Lecture Notes in Math- 
ematics, 977 (1983), Berlin- Heidelberg-New York, Springer- Verlag . 

[11] E. Plahte, T. Mcstl, W.S. Omholt, Feedback circuits, stability and multi- 
stationarity in dynamical systems, J. Biol. Sy.s. 3 (1995), 409-413. 

[12] E.H. Snoussi, Necessary conditions for multistationarity and stable period- 
icity, J. Biol. Sys. 6 (1998), 3-9. 

[13] R. Thomas, On the relation between the logical structure of systems and 
their ability to generate multiple steady states or sustained oscillations, 
Springer Series in Synergetics, 9 (1981), 180-193. 

[14] R. Thomas, Logical description, analysis, and feedback loops, in Aspects of 
Chemical Evolution, G. Nicolis ed., XVII-th Solvay Conference on Chem- 
istry (1980), John Wiley and Sons, 247-282. 

[15] R. Thomas, M. Kaufman, Multistationarity, the basis of cell differentia- 
tion and memory. I. Structural conditions of multistationarity and other 
nontrivial behaviour, Chaos 11 (2001), 170-179. 

Keywords : interaction graph, multistationarity, Jacobian matrix, global 
univalence. 



19 



